EXACT SHORT TIME DYNAMICS FOR STEEPLY REPULSIVE POTENTIALS 



James W. Dufty 
Department of Physics, University of Florida 
Gainesville, FL 32611 

Matthieu H. Ernst 
Instituut voor Theoretische Fysica, Universiteit Utrecht 
Postbus 80.195, 3580 TD Utrecht, The Netherlands 
(Dated: 1-25-04) 

Abstract 

The autocorrelation functions for the force on a particle, the velocity of a particle, and the transverse momentum 
flux are studied for the power law potential v{r) = tiajrY ( soft spheres). The latter two correlation fmictions 
characterize the Green-Kubo expressions for the self-diffusion coefficient and shear viscosity. The short time dynamics 
is calculated exactly as a function of v. The dynamics is characterized by a universal scaling function S{t), where 
T — I/ti, and t^^ is the mean time to traverse the core of the potential divided by v. In the limit of asymptotically 
large v this scaling function leads to delta function in time contributions in the correlation functions for the force 
and momentum flux. It is shown that this singular limit agrees with the special Green-Kubo representation for hard 
sphere transport coefficients. The domain of the scaling law is investigated by comparison with recent results from 
molecular dynamics simulation for this potential. 



I. INTRODUCTION 

The dynamics of fluctuations in simple classical fluids is a well-studied problem over a wide range of densities 
and temperatures. Important questions remain open at the quantitative level, but important qualitative features (e.g. 
behavior at high density, long times) have been largely resolved over the past few decades. In most cases the qualitative 
features do not depend sensitively on the form of the pair potential for interactions among the particles. An exception 
is the short time behavior of time correlation functions where the dynamics does depend sensitively on the form of the 
potential. This is due to the dominance of trajectories of pairs of particles as they traverse their common force field 
on this time scale. The most striking example of this is the difference between a continuous potential with a finite 
pair interaction time and hard spheres, for which this time is zero. The objective here is to explore this difference 
quantitatively for the case of time correlation functions characterizing the Green-Kubo transport coefficients. 

This work is motivated by the recent series of molecular dynamics studies of the same problem by Powles and 
co-workers QI^I^Q. They consider a steeply repulsive potential of the form v{r) = e{a/rY with exponent 12 < 
1/ < 1152. In the following this will be referred to as the soft sphere potential. Clearly, for asymptotically large ly this 
approaches the hard sphere potential which is infinite for r < a and zero for r > a (see Figure 1). 

It is straightforward to show that the thermodynamic and structural properties of the soft sphere fluid are contin- 
uously related to those of the hard sphere fluid, as expected. However, the corresponding relationship for dynamical 
properties is more complex. For example, the exact short time expansion of time correlation functions for the soft 
sphere fluid is a series with only even coefficients , while that for the hard sphere fluid has finite odd order contribu- 
tions as well 0. A second qualitative difference is the form of the Green-Kubo expressions for transport coefficients. 
For the soft sphere fluid these are time integrals of the flux autocorrelation function, with the flux being associated 
with some conserved density. For the hard sphere fluid it looks as if there is an additional term due to instantaneous 
momentum transport caused by configurations for pairs of particles initially at contact. Nevertheless, it is clear that 
the soft sphere and hard sphere fluids should be physically equivalent for large v. So, such apparent differences must 
be understood in an appropriate context. 

The context refers to the time scale on which the comparison is to be made. For the soft sphere potential there is 
a characteristic force range = a/v around r ~ cr, and an associated time = ^fpmajv, which is essentially the 
time it takes a particle to traverse the steep part of the potential. For large v and i > Tj^ a pair of particles initially 
separated by r « ct will have transferred an amount of momentum quantitatively approaching that for a pair of hard 
spheres. On this time scale all dynamical properties of the soft sphere and hard sphere fluid should be comparable. 
However, for t < t^, the soft sphere and hard sphere fluid are always qualitatively different in their dynamics regardless 
of how large v is taken. In effect the hard sphere fluid is always in the domain of times large compared to Ti,, as 
the collisions are instantaneous. This explains, for example, the above difference in the two short time power series 
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FIG. 1: Soft sphere potential v{r) as a function of r/a for several values ofv. 



expansions. One of the main results described here is an exact determination of the crossover behavior for the soft 
sphere fluid from t < to the hard sphere form for t > Tu for large v. 

The detailed analysis given here is made possible by the simplifications that occur for repulsive power law potentials 
at large v. This entails both a limitation to times t k and restricted spatial domains a — < r < a + r^. The 
relevant other time and space scales are the mean free time Ie and the mean free path Ie- Although both of these can 
be quite small at high densities (e.g., Ie < c) they are insensitive to v. This means that for sufficiently large v there 
is a separation of both time and space scales, << Ie and « Ie- As a consequence the dynamics to be studied 
reduces to pair dynamics since force ranges of different pairs will not overlap and times are much shorter than that 
for sequences of pair collisions. The problem of evaluating the pair dynamics on this time scale has been solved for 
the velocity autocorrelation function (VACF) for three-dimensional soft spheres by de Schepper His analysis can 
be extended to other time correlation functions and different dimensions as well. It is found that the time domain 
for t «tE and «tE is described by a universal scaling function S{t) depending on time only through t = t/Tn 
and otherwise independent of density, temperature, and potential parameters. As it turns out, the time dependence 
of all Green Kubo correlation functions are related to a single scaling functions S{t). The domain of validity for this 
scaling law is studied by comparison with the three-dimensional simulation results of Powles et al. at v = 1152 and 
at a packing fraction of ^ = 0.3. A more complete comparison at v = 1152, at higher densities, and for more general 
time correlation functions will be given elsewhere 

The initial values of correlation functions of fluxes involving the force are proportional to ly for large v • The 

combination of vS{t/T,j) becomes proportional to a Dirac delta function in t in the limit ^ cx3 at flxed t. This is the 
domain t > t^, for which the hard sphere limit is expected. To explore the crossover to hard sphere behavior and the 
implications of this delta function a detailed description of the Green-Kubo relations and time correlation functions 
for the hard sphere fluid are included here as well. The dynamics for hard spheres is no longer described by forces. 
Instead there are straight line trajectories for all particles until a pair is in contact. Instantaneously, the pair exchanges 
momentum according to an elastic collision and proceeds along the new straight line trajectory. The generators for 
this dynamics (Liouville operators) involve binary collision operators rather than forces jlJl- These differences from 
the dynamics for continuous potentials lead to the qualitative differences in the Green-Kubo relations and associated 
correlation functions mentioned above. Generically, the time correlation functions C{t) approach in the limit as v —^ oo 
to C{t) = £oo(5+(t) -I- 6C{t). The first term represents a singular part, which is delta-correlated in time, and a regular 
smooth function 6C{t). Here S^{t) is the delta function normalized to unity over the positive time axis. The coefficient, 
Coo, vanishes for any continuous potential but is non-zero for hard spheres. Consequently, the Green-Kubo formulas 
for the transport coefficients, £, in the hard sphere limit have the form, C = Coo + fimt_»oo fimy^oo /g ds6C{s). 

This limiting behavior is confirmed analytically below and shown to be consistent with the simulation data. Al- 
though the literature on hard sphere fluids is large, the detailed forms for hard sphere Green-Kubo relations and 
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discussion of these differences do not seem to have been given before. Instead, an equivalent form of Helfand relations 
has been used in the literature both for theoretical analysis 1111 1141] , and for computer simulations ll5|| , since this 
form does not involve the forces explicitly, and is valid for both the soft sphere and the hard sphere fluid. Here the 
Helfand relations are taken as the starting point for derivation of the hard sphere Green-Kubo relations. 

A preliminary report on the issues considered here has been given by one of us (16| , restricted to the case of the shear 
viscosity and without calculating the scaling function. There is a substantial literature on the dynamics for continuous 
potentials on time scales comparable to time , typically associated with memory function models |l7l | . This domain 
is important for many conditions in neutron scattering, spectroscopy, and short-pulse laser experiments where high 
frequency domains can be accessed. The attention here is more focused on the limiting form of the dynamics for the 
soft sphere potential and the associated universal scaling properties. For a given potential, e.g. Lennard- Jones, the 
short time dependence will be potential specific and the results obtained here have little or no relationship to more 
general potential forms. 

The analysis here is exact in the limits considered. The objective here is to clarify the differences mentioned above 
for soft sphere and hard sphere interactions regarding the dynamics of the associated correlation functions. To allow 
this broad scope of discussion, attention is limited to the auto correlation functions for the force, which exposes 
most directly the scaling function at short times, and those for the shear viscosity and self-diffusion coefficient. A 
similar analysis of the thermal conductivity and bulk viscosity has been carried out and will be present separately 
The plan of the paper is as follows. In the next section the Green-Kubo relations for the shear viscosity and self- 
diffusion coefhcient are recalled for the case of smooth (diffcrentiable) potentials. The fluxes defining the associated 
correlation functions are identified as the sum of kinetic part (k) and a coUisional transfer or potential part (v), 
and the correlation functions are decomposed into the corresponding contributions from these components. Also in 
this section the equivalent Einstein-Helfand formulas are recalled. These formulas are applicable to both smooth 
interactions and hard spheres. The generators for hard sphere dynamics are given in Appendix A and applied in the 
Helfand formulas to derive the Green-Kubo formulas for hard spheres discussed in Section HI. For later comparison 
with the limiting forms for the soft sphere fluid, the short time behavior of the hard sphere correlation functions is 
calculated in Section IV. The short time dynamics of the soft sphere fluid is addressed in Section V, where the scaling 
properties of the velocity autocorrelation function and force autocorrelation function arc discussed. The dominant 
short time contribution to the autocorrelation function for the shear viscosity, which is the stress autocorrelation 
function, is given by the same scaling function. These results are then compared with the MD simulation results of 
Powles et. al. in Section VI and final comments are offered in the last section. 



Our goal is to study the short time behavior of time correlation functions C{t) = {P/V){J{0)J{t))^ of microscopic 
fluxes J that enter in the Green-Kubo formulas for the transport coefficients in classical fluids with smooth interactions. 
These fluxes contain in general a kinetic part, J*^, and a potential part, J", i.e. J = -\- J". In this paper we only 
illustrate the general results by considering the most simple cases, being the velocity autocorrelation function (VACF), 
the force autocorrelation function (FACF), and the stress autocorrelation function (SACF). The time integrals over 
the VACF and the SACF determine the self-diffusion coefficient and shear viscosity, respectively. 

For smooth interactions the time correlation functions are regular at the origin, i.e. they can be expanded in powers 
of . In the limit where the repulsive part of the potential becomes very steep, and approaches the hard sphere limit, 
singularities develop at i = 0. The type of singularities is quite different for the different afc— parts of the correlation 
functions, C^t) = {(3/V){J''{0)j''{t)) ^ with ab ^ {kk,kv,vk,vv}. 

We start with the Green-Kubo formula for shear viscosity expressed in the grand canonical ensemble, characterized 
by a temperature T = l/ksfS, a chemical potential, and a volume V. Here averages are denoted by (• • • )q. Then, 



II. TIME CORRELATION FUNCTIONS FOR SMOOTH INTERACTIONS 
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where r]vit) is a time integral over the time correlation function, 

c,{t) = ip/v)(mj{t))^ = cf{t) + 2cf (t) + a™(t). 

Here C^"" {t) = Cli'^ (t) are equal, and the fluxes are given by 
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Similarly, the Green-Kubo formula for the self diffusion coefficient D is given by 

D — lim lim Dv{t) = lim lim / dsCuis), 
where C]j{t) is the autocorrrelation function for the velocity of one of the particles, 

A self-correlation in the grand canonical ensemble has to be understood as {ai{0)ai{t)) ^ = {J2i^ii^)'^ii^))o/{'^)o- 
Finally, we will also consider here the autocorrelation function for the total force acting on a single particle 

Cpit) = {(3/m){F^.,{0)F,,{t))^. (II.6) 

where Fix — J2i>i Pii,x is the force on particle 1. The FACF is simply related to the second derivative of the VACF. 
Its time integral does not give any transport coefhcient, but vanishes in fact. However, it is included in this study since 
it has a time dependence closely related to that of the SACF. For simplicity it will be referred to as a "Green-Kubo 
correlation function" along with Coit) and C^(t). 

Next we consider the Einstein- Helfand Formulas. These formulas 13, Tsj for the transport coefficients are the 
analogs of the Einstein formula for the self diffusion coefficient D in terms of the second moment of the displacements, 



(11.4) 



Dv{t) = hiii^it) - ^m')o = /o ds{vix{0)vix{s))^. {11.7) 

The first equality is the Einstein- Helfand form while the second equality is the Green-Kubo form. The latter follow 
directly from the identities x = Vj. and x{t)~x{Q) — dsx{s), and the stationarity of the VACF. The Helfand formula 
for the shear viscosity is given by an analogous moment formula, 

M = mVixTiy, (II. 8) 

i 

where M is a Helfand moment. The equivalence of ljll.8|) with the Green-Kubo formula can be established along the 
same lines as in (III.7|) by observing that the flux is given by, 

J = M = {Af, H} = LM, (II.9) 

where H ^ K + V is the Hamiltonian, and {M, H} is a Poisson bracket. The last equality defines the Liouville 
operator for smooth potentials. Expressions in terms of a Liouville operator are of interest as we want to consider also 
interactions between hard spheres. Such interactions can not be described by a Hamiltonian being a smooth function 
of the relative distances between the particles. However, there exist in the literature pseudo-Liouville operators, that 
generate the hard sphere dynamics inside statistical averages 

m The Helfand formulas are given in terms of the 
moments M which do not involve the force that becomes singular for the hard sphere limit, in contrast to the fluxes 
in the Green-Kubo representation. As indicated in the next section, use of the hard sphere pseudo-Liouville operators 
in ljll.91) provides the correct definition for the microscopic hard sphere fluxes and the corresponding hard sphere 
Green-Kubo representation using (|II.9p . 



III. TIME CORRELATION FUNCTIONS FOR HARD SPHERES FLUIDS 



A. Einstein-Helfand and Green-Kubo Formulas 



The corresponding time correlation functions for hard sphere fluids, denoted by C^''{t), are singular at t — 0. They 
behave quite differently at short times from those for smooth interactions, which are regular at i = 0. For instance, 
the VACF for hard spheres decays exponentially on the time scale of the Enskog mean free time ts, and it is even in 
t. So, it is singular at i = with a jump in the first derivative, and a delta function in the second derivative. 

In order to investigate what happens in the limit of very steep repulsion, we want to first calculate the corresponding 
hard sphere results, which have only been studied in the literature for the VACF and the incoherent scattering function 
dill. The problem is that the fluxes J in ljll.3|l contain the forces Fy , which are ill-deflned for hard spheres, and 
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there is no obvious way to extend the usual Green-Kubo formulas to hard sphere fluids. However, as noted above, the 
equivalent Einstein- Hclfand formula (|11.8|l involve only momenta and energies which remain well-defined in the hard 
sphere limit. 

The proper way to formulate hard sphere dynamics for an equilibrium time correlation function {A{0)B{t))^ is 
to introduce the forward (+) and backward (-) generators e**^± with t > 0. They generate the trajectories in T— 
space, A{t) — e^^^^ A{0), for an arbitrary phase function inside a statistical average. The generators involve the 
pseudo-Liouville operators L± and the binary collision operators T± defined in Appendix A (see Ref [l3j). where we 
also explain how to express the Einstein-Helfand correlation function for rj in terms of these hard sphere generators, 
and subsequently into a Green-Kubo formula. The result is, 

wit) = f dsC^is) - 7?oo + / dsSCrjis), (111.1) 
Jo Jo 

where the limits in still have to be taken. Here we have introduced, 

a,(i) = VocS+it)+dCr,ity, SCriit) - (/3/V^)(J_e*^+J+)„ 
J± ^ J*' + Jl= L±M; J*" = LqM = ^ mvi^v^y 

i 

J± = T±(ij)M = ± ^ T±{ij)\mgij^^rij,y 

i<j i<j 

Voo = -iP/V)(ML+M)^. (111.2) 



The delta- function is normalized as dtS^{t) — 1. The kinetic part of the flux in ljni.2|) is identical to that for 
smooth interactions (|11.3|l . For the time being the coUisional part, involving T— operators, will be kept in the schematic 
form above. 

Next we consider the instantaneous viscosity 7700, which is defined as an equilibrium average, and vanishes for 
smooth interactions ( i.e., when L-|_ is replaced by L). For the hard sphere fiuid 7700 can be expressed in terms of the 
hard sphere pair distribution function by combining the expressions for 7700 and L+M — J+ in (|III.2|I . where we have 
changed to center of mass, {Ry, Gy — ^(vi -I- vj)}, and relative phase variables, {ry, gy = Vi — vj}. The result is 

r/00 = -\l3{mn)-'x! dv{{g^Ty T+(12) g,r,)), (HLS) 

where ((• • • )) denotes a Maxwellian velocity average over all particles involved. Moreover, the r— integration can be 
carried out because the operator T+ contains a factor 5'^'^\r — aa). Consequently x = the hard sphere pair 

correlation function at contact. The remaining integrals are d— dimensional generalizations of the collision integrals 
as appearing in the Enskog theory for hard sphere fluids (See Chapter 16.8 of Ref. Performing the a— and 

velocity integrations yields finally, 

mncr^ d 



The Doric pi w, defined in ljni.4l) . has been chosen such that it reduces for d = 3 to the same symbol as used in the 
classical Enskog theory, as presented in Chapters 16.51, 16.52 and 16.6. of Ref. [l^. In the above formula the Enskog 
mean free time t^; is given by, 

tE = y/TTt„/2dhnx = ^tal2dl\ (HI.S) 

where — \J (5raa and h is the excluded volume, equal to half the volume of the d-dimensional interaction sphere of 
radius a. 

There are two conspicuous differences between the formulas for hard spheres and for smooth interactions. First, the 
pseudo-fluxes J± in the time correlation functions are different depending on their position relative to the generator 
6*''^+. We also note that J+e"*^- J_ is an equivalent order. Second, there is the instantaneous contributions, 7700, 
which is vanishing for smooth interactions. 

Before closing this section we point out that the hard sphere time correlation formula C,j{t) in ljni.2p can also be 
split into C^''{t) with (ab) — {kk, kv, vv} by splitting the flux as J± — + (see (jHI.2|l '). The form of the kinetic 
part J'^ is identical to that for smooth interactions. The remaining collisional transfer part is different. 
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B. The hard sphere FACF and VACF 

In the same way as described above for C,,(s), the Helfand representation for the force autocorrelation function can 
be used to obtain the equivalent hard sphere form 



TO 

7oo = -(3m{vixL+vix)^::^ iA/V^U = 2/dtE. (III.6) 

One directly recognizes the expression for 700 as the opposite of the initial slope of the hard sphere VACF, which 
has been calculated exactly along the same lines as (|IIL4|) (see Ref.0|). We also note that 700 vanishes for smooth 
interactions. 

Finally, the VACF for hard spheres is much simpler to obtain since there is no explicit dependence on the forces. 
Consequently, only the generator for the dynamics has to be changed, leading to 

Coit) = (t;i.(0)e*^+t;i,)„ (III.7) 

This completes our identification of the hard sphere Green-Kubo time correlation functions (7^ (t), Cp (t) , and Cd (t) ■ 



C. A new relation for the hard sphere fluid 

An interesting consequence for the hard sphere FACF follows by explicitly integrating (|II.6|I over t using Newton's 
law to find for, say, soft spheres, 

lim / dsCpis) = f3 lim (^^,(0)(«i,(i) - «i.(0)))„ = 0, (IIL8) 

t — >oo Jq t — >oo 

because positions and velocities become uncorrelated in this limit. Taking the limit of the exponent of the power law 
interaction ^ 00 suggests that this relation also holds in the case of the hard sphere correlation function C{t) in 
(inra . i.e. dsC^is) = 0. As a consequence we obtain the following relation, 

7oo = -- / dt(Fi,_e*^+Fi,+)„. (III.9) 

The left side of the last equation has been calculated exactly above. The right side is a complex dynamical quantity 
involving the entire time evolution of the system. Thus we have obtained a rare "zero frequency" sum rule for the 
hard sphere fluid. It should be noted however that equating ljIII.8|) with the corresponding integral over C, implies 
an interchange of limits, i.e. limt^oo and lim,y_»oo, which is presumably allowed. However, interchanging the limits, 
limt_,o and lmi^_,^, is not allowed, and the consequences of this non-uniformity is in fact the main subject of this 
paper. The full implications of the new relation ljIII.9|) are not clear at this point. 



IV. EVALUATION OF HARD SPHERE PROPERTIES 

Having obtained Green-Kubo correlation functions C°^(<) for smooth interactions and C^^{t) for hard sphere in- 
teractions, we now discuss the structure of the short time behavior of these functions. After a brief introduction we 
consider in this section the hard sphere case, and in Section V the case of smooth interactions, in particular soft 
spheres. 

First, some perspective on the differences between these two cases is given by listing in Table I the qualitative 
structure of the short time behavior for smooth and hard sphere interactions. We first note that all Green-Kubo-type 
correlation functions are even functions of time. If they are regular at t = 0, as is the case for smooth interactions, 
then they can be expanded in powers of t^. For hard sphere systems the correlation functions are also even functions 
of t, but they are singular at the origin. 
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Short time behavior of time correlation functions 
for smooth interactions and for hard sphere interactions. 



Mi 


Cf{t) FOR SI 


C^^''{t) FOR HARD SPHERE 


ab 


kk 


kv 


vv 


kk 


kv 


vv 


D 


0(l) + 0(t2) 






0(l) + 0(|i|) 






F 






0(1) 






S{t) 


V 


0(l) + 0(t2) 




0(1) 


0(l) + 0(|i|) 


0(1) 


m 



Table I shows the schematic structure of the selected correlation functions for smooth interactions and hard sphere 
fluids, where the terms 0(i") and 0(t") with n = 0, 1, 2 are non-vanishing terms of order for small t. The results 
for smooth interactions are well-known in the literature. Regarding hard spheres we first note that at small t only the 
leading order terms (initial values) in C^'' (t) and C^'' (t) are equal. The entries in the remaining columns for smooth 
interactions and hard spheres are all different. The hard sphere results are obtained from l|III.2|l for Cjj{t), (jIII.6|) for 
Cpit), and ljIII.7|l for Coit). All contributions involving a single T— operator are non- vanishing, as in the hard sphere 
entries on location (/i, kk) and (/i, kv). The hard sphere entries on location (/x, vv) involve two T— operators, and are 
more complicated. They will be discussed later. 

Inspection of Table I shows that the results for smooth interactions and hard spheres are indeed very different. 
The goal of this section and the next is to calculate the entries in the table, and study in a quantitative manner 
the crossover of the correlation functions to hard sphere interactions from smooth, but steeply repulsive power law 
interactions, v{r) ~ 1/r'^. This will be done for times t, short compared to the Enskog mean free time Ie- 

In the remaining part of this section the short time behavior of the correlation functions for the hard sphere fluid 
in Table I will be calculated, i.e. the initial values and initial slopes. We start with the (kk)-correlations, and include 
the VACF Coit) as the most typical one. In the sequel we restrict ourselves exclusively to t > 0, to avoid possible 
confusion regarding the deflnitions of the T— operators. 

Coit) = {vi4l+ tL+ ■■ -jvi^) ^= {l/(3m){l--/Dt + ---} 

Cf{t) = (j'={l+iL+...}J^-)„^(n//?){l-V + ...}. (IV.l) 

As explained in the previous subsection, we restrict ourselves in the small-t expansion to terms that are at most linear 
in the T— operator. Terms of 0{{tT)'^) have been neglected in ljIV.l|) . Initial values and slopes can be evaluated, and 
yield 

ID = Joo=-{VxL+V^)^/{vl)^ = 2/dtE 

7, = -{j'L+j')J{{j')')^^i/{d + 2)tE (IV.2) 

with J'= defined in (p3|) . _ 

Next we consider the (kv)-cross-correlations C'^'"(t) = {/S /V){J'' exp[tL+]Jl) ^. As Jl = J+- j'' in l|lIL2)l itseff is 
already linear in T, we can only calculate its initial value exactly. 

Of (0 = if3/V){JHl + ■■■} E^<, T+{ij)\mg,,^xn,,y), 

= i/3(mn)2x/dr((5.g,r+(12)g,r,)) = (l)^ . (IV.3) 

Finally we discuss the (vv)- or coUisional transfer correlations for hard spheres in Table I. It is instructive to first 
compare 0™(i) for smooth power law interactions with C^^{t) for hard spheres. The correlation function C^'"{t) for 
smooth interactions develops in the hard sphere limit a strong delta function - type singularity, as represented by the 
first line in (|III.2(I - and similarly in l|III.6p for Cpit). The remaining part, 6C^'"{t), containing the pseudo-fluxes J^., 
represents in fact the regular part, that approaches a finite limit as t — > 0. Indeed, the short time behavior of SCpit), 

or equivalently Coit), has already been analyzed in great detail in the literature (see Ref.0|). There it has been 
shown that the hard sphere correlation functions related to SCoit) or SCpit), and containing two T— operators, are 
smooth functions of time near t = 0, which indeed approach a finite non-vanishing value. The explicit evaluation of 
these contributions is much more complex than performing simple binary collision integrals, as we have been doing 
in the previous sections. The reason is that these contributions are coming from the overla ppi ng part of uncorrelated 
binary collisions (12)(13), and from renormalized ring collisions of the form (12)(13)(23) [lal, which involve in fact 
three instead of two hard spheres. However, for the purpose of this paper the fact that these values are finite, is 
sufficient. 
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V. SCALING FORMS IN SOFT SPHERE FLUIDS 
A. Force autocorrelation function (FACF) 

In the previous section we have calculated the leading short time hard sphere results. These are needed to compare 
and identify the limiting results for soft sphere correlation functions on time scales that describe the crossover to 
hard sphere behavior, i.e. the crossover from the initial time scale, where the detailed shape of the interparticle 
interaction matters, to the kinetic time scale tE, where only asymptotic scattering properties matter. 

This will be done for soft spheres, represented by the repulsive power law potential v{r) = t(a jrY ■ In this case the 
FACF and SACF turn out to be proportional to a scaling function SitjTy), where Ty — t^jv is the mean time that a 
particle needs to traverse the steep part of the potential, and ta — \J fima the mean time to traverse the total hard 
core diameter a . These time scales are only well-defined and well-separated for large values of exponent v. At high 
densities there is another relevant time scale, the mean free time between collisions, tE, which can be estimated for 
sufficiently steep repulsion by the Enskog mean free time, tE, which is proportional to t„ /hnx- Given the large values 
of the two- and three-dimensional pair correlation function x ^-t contact, tE and t^j can be of the same magnitude (for 
example in two- and three-dimensional hard sphere systems at packing fractions around 30%), ov tE may even be an 
order of magnitude smaller than t„ (for example in typical neutron scattering experiments on liquid Argon). Hence, to 
describe the crossover from the initial Ti,— scale to the kinetic scale tE, the initial scale must satisfy <C minji^, t„}. 
These estimates suggest that the exponent v should be rather large at high densities. The example of the next section 
with packing fraction 0.3 and v = 1152 satisfies these constraints, as will be illustrated in Section VII. 

Our analysis begins with the FACF in a d-dimensional system, defined in (|II.6|I . It is an even function of time t 
and regular at the origin for finite values of the exponent v , i.e. it can be expanded in powers of t^. However for 
very large v, the initial value Cf(0) — > oo, implying that the function is singular at the origin. It is the goal of this 
subsection to analyze the dominant small-i singularity of Cpit) in the hard sphere limit and to describe the crossover 
of the FACF from soft to hard sphere behavior. 

A study of this short time crossover problem can be carried out following the work of de Schepper |^ for the 
VACF. His analysis is based on a perturbation expansion of e*^ in powers of (tL)^, where the Liouville operator 
contains the force Fjj. He has shown that the most dominant contributions for large u are obtained by keeping in 
each order in (ti)'^ only terms involving forces between a single pair {ij), and finally resumming these contributions. 
Here we exploit this result and calculate directly the entire pair contribution. The basic physical idea is that the 
autocorrelation function (Fi2(0)Fi2(t))Q of the pair force F12 = — Vw(r) ~ 0{i') controls the short time dynamics on 
the time scale t^, and the time evolution is controlled by two-particle dynamics since t^, <C tE- More explicitly, the 
dominant short time contribution to CF{t) for 00 is, 

CF{t) =i {l3n/mV) j dRdrdgdG (f>{g)(f>{G)g'-^'> {r)F^{r)e*^^^ F^{r) 

= (/3nxM)(^i2,xe*^'-Fi2,,)''^. (V.l) 

The second line defines a two-particle average over positions and velocities, the latter one with Maxwellian weights. 
A change of variables from {ri, Vi, r2, V2} to relative and center of mass coordinates {r, g, R, G} has been made with 
the replacement L12 = Lr + G ■ Vr, and n — {N)^/V. For large v the r— integrand is sharply peaked around r = a. 
Consequently the pair distribution function g^^'^ir) can be replaced by its value at contact, x = ff^^"* ("'+)■ 

The analysis is now reduced to a one body problem in the soft sphere potential. The detailed calculations are still 
rather technical, and will be published elsewhere. The final result is 

Cpit) ~ iloo /T,)S{t/T,). (V.2) 

where the crossover function S{t) is found as, 

/ J \ 2 /.oo 

5(r) = ^— j 2t J dye-y\^cothTy. (V.3) 

Interestingly, the crossover function is independent of the dimensionality d, and depends on the exponent only 
through the scaling variable r = vt/tc- All explicit j/— dependence in (jV.2p is accounted for in the overall factor 
^'contained in I/t^. 

The crossover function S{t) can be expanded, both at short and at long times, in a convergent infinite series with 



FIG. 2: Figure (a) shows the exact scaling function 5(r) (solid line) together with the phenomenological one l/cosh(r^) 
(triangles), as discussed in the text. The insert shows the same on a logarithmic scale, in order to visualize the differences 
between the exact and phenomenological one at large r. Figure (b) shows the scaling function J-{r) (solid line) for the VACF 
together with its large— r asymptote (dashed line). 



known coefficients, of which we quote the leading terms, 

(^<'^) (V.4) 
^' \^T-5 + 0(r-^) (r>7r) ^ ' 

A numerical evaluation is shown in Figure 2. 

The scaling property was first noticed in the MD simulation results of Powles and collaborators Q, |^ |^ Q . In that 
work a phenomenological crossover function cx 1/ cosh(T'\/2) has been used, which decays exponentially. The insert 
in Figure 2 also shows that it is correct at short times, but decays too fast at longer times. 

The short time behavior for the soft sphere fluid (fixed, but large v) occurs for r = I/ti, « it and therefore is 
obtained from the small— t behavior in 



CF{t) = ^ [y^ [1 - [t/T^f + 0{{t/T,t)\ . (V.5) 



To describe the crossover to the hard sphere fluid we observe that 

lim —s(—]^ 6+{t) (V.6) 
can be considered as a delta function on the interval (0, oo), i.e. 

/•OO /'OO 

/ dt6+{t)= lim / dT5(T) = l. (V.7) 
Jo ''^°°Jo 

Its value at f = is infinite, and vanishes at i 7^ 0, because 5(0) and S{oo) = 0. The hard sphere behavior 
{ly 00) at small, but fixed t , is obtained from the large— t behavior of S{t) in (jV.4|l and l|V.6|) . and yields for 

t/Ti, > TT or ly > irta/t, 

Cpit) = ^ [y^ [s+iy) + \l^\T^/tf + 0{{t/T,y)^ . (V.8) 



The first term inside the curly brackets represents the dominant contribution to the hard sphere FACF in ljIII.6|) at 
short times, and the second term represents the dominant correction to the limiting result. The identification of the 
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crossover function and demonstration that it generates a delta function singularity is one of the main results of this 
work. It provides the basis for understanding the connection between the soft sphere VACF and the singular form 
of the hard sphere VACF shown in l|III.6(l . It is noted that the coefficient of the delta function arising from (|V.2|I is 
exactly the same as that of (jIII.6ll coming from transformation of the Helfand formula. 



B. Stress Autocorrelation Function (SACF) 

Application of this analysis to the collisional transfer or vv-parts of the correlation functions for the shear viscosity, 
bulk viscosity, and thermal conductivity leads to exact results for the short time crossover function of the form ljV.2p 
with a different pro-factor, but with the same scaling form S{t). In the hard sphere limit this scaling form approaches 
according to (|V.6|I again to a Dirac delta function. The kk- and kv-parts of the time correlation functions at short 
times are less singular than the vv-parts. So the dominant short time singularity of the full time correlation functions 
is contained in the collisional transfer (vv) terms. We only quote explicitly the result for the stress-stress correlation 
function, 

C,{t) c C:;{t) c (r?oo/T.)5(Vr.). (V.9) 

Again, the delta function singularity and its prefactor associated with this result agree exactly with the singular part 
of JTTO). 



C. Velocity Autocorrelation Function (VACF) 



Finally we consider the VACF Cnit) for the soft sphere fluid, as defined in (|II.5|I . It is regular at < = with 
Coit) = (1//3to){1 + 0{t'^)} according to Table I. On the other hand, for the hard sphere fluid Coit) = (l//3ra){l - 
7oo|i|/io- + •••} is linear in t at small positive t. Again, computing the dominant contribution of the pair force 
autocorrelation function with pair dynamics we recover the result of , 

Coit) = {1 + 7ooT-..F(r)} = Cz5(0) (l + ^^(r)l (V.IO) 
with A = bnx and r = t/r^, and the new crossover function is 

/•oc 

T(t) = 2 dye-y\^{l-TycothTy). (V.ll) 



The function J-{t) is simple related to S{t) by 

{S/dT'')T{T) = -S{t). (V.12) 

This is consistent with the exact relationship of the VACF and the FACF, Cpit) ~ — PmC oit) ■ The initial conditions 
on T{t) as t can be read off from location (D, kk) in Table I to be — and ^'(0) = 0, which determines 
the two integration constants. 

The small and large r— expansions of J-{t) are for positive r, 

^ ^ - I -r + i - i^r- + 0{r-^) (r > '^''^ 

The short time behavior in the soft sphere system at fixed, but large v is now found to be, 

Coit) =^ {1 - Ht/rf + 0{{t/T,Y)] . (V.14) 

This is the explicit form of the short time behavior of Cnit), listed in schematic form in Table I on location {D, kk) 
for the smooth interactions case. Similarly we find for large t/T^, behavior (with t fixed, but small), 

Coit) ^ {1 - (4A/ V^)[t/r. - + 0{r^,)} , (V.15) 

where the first two terms inside the curly brackets represent the behavior of the VACF for hard spheres at small t, 
as given schematically in Table I, and explicitly in ljIV.l|) . Note that the dominant correction of 0{i'~'^) to the hard 
sphere result is independent of t, and decreases very slowly with v. 
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FIG. 3: Force autocorrelation function (FACF) as a function of t = I/ti, (solid line) for packing fraction 0.3 and v = 1152, 
compared with MD simulations (circles). The insert shows the same except on a logarithmic scale, in order to make manifest 
that the asymptotic theory breaks down at large t. 



VI. COMPARISON WITH SIMULATIONS 



During the past five years Powles and collaborators have carried out a systematic study of the short time dynamics 
for the soft sphere fluid by molecular dynamics simulation. They report resuhs at two packing fractions 0.3 and 0.45, 
corresponding to a moderately dense fluid, for the range of exponents v = 12,24,36,72,144,526, and 1152. This 
is done for the correlation functions occurring in the Grccn-Kubo expressions for shear and bulk viscosity, thermal 
conductivity, self-diffusion coefficient, and the force autocorrelation function. These results provide extensive data to 
test and interpret the theoretical predictions made here. As only partial results have been given here, for the FACF, 
SACF and VACF, the comparison will be limited to these cases. Furthermore, attention will be restricted to the 
packing fraction 0.3 and i/ = 1152. These arc the best conditions for the separation of time scales <C minjts, fo-}- 

We first observe that our theoretical results are asymptotically exact, i.e. ly has to be large enough such that 
on the one hand the gap between and the minimum of {is, to-} is sufficiently large to test the predictions for 
asymptotically large r = t/Ti^, and on the other hand t has to be short enough, i.e. t << ts, such that the dynamics 
of the isolated pairs describes the full time evolution of the TV— particle system. An estimate at ^ = 0.3 and = 1152 
gives: ts — 125ri, ~ 0.099io-, illustrating that the time scales are well separated. 

It is useful to anticipate the results of the comparison between the short time scaling behavior determined here 
and that from the simulations. Since the results here are asymptotically exact, there must be agreement in all cases 
where Ti, « {tE,tcr} and t « ts- There will then be the crossover domain where the short time form begins to 
fail. However, the long time form of the scaling functions represents the "short time" dynamics of a hard sphere fluid 
(fixed small t, large v), as indicated in the previous section. This is the domain for which the VACF exhibits behavior, 
linear in t. 

What happens for larger t, where the decay of the VACF becomes exponential, rather than linear (say for t > OMe), 
is strictly outside the predictions of the present theory, and no exact results are known for these larger time intervals. 
Of course the results may still be compared with the Enskog theory for hard sphere fluids, but wc do not consider 
approximate results in the present article. Similarly, the singular part of the FACF and SACF will be very small on 
this kinetic time scale, and no predictions on the regular parts of these functions, i5Cf(t) and SCn{t), can be given in 
the context of the present theory. We note that these regular parts contain, among other, the full kk- and kv-parts of 

c,{t). 

Consider first the FACF. Figure 3 shows that the agreement between the short time crossover function and MD 
simulation is excellent at very short times. The insert shows the same data on a logarithmic scale to discover that 
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FIG. 4: Stress autocorrelation function (SACF) as a function of r = f/r^ (solid line) for packing fraction 0.3 and i/ = 1152, 
compared to MD simulations (solid circles). The insert is a logarithmic plot of the same data, to show that the asymptotic 
theory breaks down for large r values. Note that the MD data (open circles) for the coUisional (vv) part of the SACF agree 
with the theory over a larger time interval than the full SACF, as is to be expected. 

deviations occur for t > 3. Figure 4 shows a similar representation for the SACF. Also shown in the insert for this 
figure is the potential part of the correlation function. It is seen that the agreement between the crossover function 
and the potential part is again excellent up to about t = 4. 

The VACF shows a smoother crossover to the hard sphere behavior. Figure 5 shows the good agreement for 
T < 4, and also shows clearly the crossover to the hard sphere form, as follows from the long time expansion. The 
extrapolation of this hard sphere form to t — gives a value different from 1, due to the crossover to soft sphere 
behavior near t = 0. As v increases this short time domain goes to zero and the offset for the hard sphere form goes 
to zero, as shown by the lower dashed curve. The large— r asymptote, y = t — shows a horizontal shift ^^/tTj 

as follows from the asymptotic expansion (jV.15|l . 



The exact short time dynamics has been calculated for the force, stress, and velocity autocorrelation functions for 
the special case of a power law potential with large negative exponent. The analysis applies for times much shorter 
than the mean free time between pair collisions ts, but includes times less than and greater than the time t^, = tu/v, 
which is the time a particle needs to traverse the steep part of the potential. This latter interval around t^, is referred 
to here as the crossover domain because the dynamics is characteristic of fluids with smooth interactions at the shorter 
times {t < T^), but becomes more representative of a hard sphere fluid at longer times {t^ « t « Ie)- The results 
can be all expressed in a single universal scaling functions, S{yt/t„), which is the autocorrelation function of the pair 
force, (Fi2,2;(0)Fi2.x(0)o times t « tE- It is independent of the dimensionality of the system. 

Thus, identification of the crossover dynamics allows several conceptual and quantitative questions to be addressed 
regarding the relationship of fiuids with steeply repulsive interactions to those of the hard sphere fiuid. To pose such 
questions it is important to have detailed results for the hard sphere fluid. There are complications due to the singular 
nature of the interactions. For example, the correlation functions associated with the Green-Kubo expressions for 
transport coefficients involve fluxes with forces that are ill-defined for hard spheres. Similarly, the usual Hamiltonian 
dynamics for these correlation functions is no longer applicable for hard spheres. 

The first part of our analysis here was to show how these problems could be handled, with well defined expressions 
for the hard sphere fluid. First, the Helfand representation for transport coefficients was recalled as an equivalent 
alternative to the Green-Kubo forms. It was shown in Appendix A how this representation could be used, together 



VII. CONCLUSIONS 
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FIG. 5: Velocity autocorrelation function (VACF) as a function of r = t/r„ for packing fraction 0.3 and v = 1152 (solid circles), 
compared to MD simulations (solid circles). Also shown is the same asymptote as in figure 2(b), parallel to the exact initial 
slope 7_D and shifted by the constant "offset" in Eg. lV.lSt . 

with an appropriate pseudo-Liouville dynamics to arrive at the corresponding Green-Kubo expressions for the hard 
sphere fluid. Interestingly, these differ in form from those for smooth interactions by constant contributions due to 
instantaneous coUisional transfer contributions from initial configurations of particles at contact. This implies that 
there appears in the hard sphere limit (y oo) a delta function singularity in the correlation functions starting 
with smooth interactions. The exact short time behavior of the velocity autocorrelation function and kinetic-kinetic 
parts of the Green-Kubo time correlations functions, are described by the scaling function !F{vt/t„), which satisfies 
<P'T{j)ldT'^ = —S(t). In the hard sphere limit these functions have cusp singularities, with exact short time expan- 
sions, depending linearly on the magnitude of i, while the corresponding functions for the smooth interactions are 
analytic in time. 

A primary result reported here is the identification of a single crossover function S(t) that describes the short time 
dynamics of all of these time correlation functions for smooth interactions. The above anomalous differences between 
smooth and hard interactions can be understood qualitatively and quantitatively from the properties of this function. 
For the stress autocorrelation function it yields a delta function singularity showing the consistency between the 
different Green-Kubo forms for smooth and hard interactions. It also shows the non-uniform nature of the short time 
limit with respect to the steepness of the potential, explaining the differences between the analytic and non-analytic 
forms of the short time expansions for smooth and hard interactions, respectively. 
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Appendix A: Pseudo-Liouville and Binary Collision Operators 

The time evolution of any phase function A{t)^ evolving under hard sphere dynamics with diameter tr can be 
described for positive times (i > 0) by the forward generator = e*-'"+A(0) with the pseudo-Liouville operator 
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The equilibrium time correlation functions necessarily satisfies the relation {A{0)B{t))^ = {B{0)A{—t)) ^ because of 
stationarity. This requires in the hard sphere case a backward generator e~*^"(i > 0), where 

L± = Lo±E.<,r±(*j) 

T±iij) = a''-' f^^da |g,, • a\ 6{r,, - <j)% - 1) . (A.l) 

The superscript on the ct— integral denotes the constraints (Tgy ■ > 0). The 6cr— operator is a substitution operator, 
executing the collision laws, i.e. 

ha-Vj = v^- = Vj + (gy • a)a. (A.2) 
The stationarity of the time correlation function under hard sphere dynamics for positive t is expressed as, 

{Be^^+A)^ = J dTpoBe*^+A^ j dV poAe'^^- B = {Ae-^^- B) ^, (A.3) 

where the operators L± and e**^± inside the thermal averages (• ■ ■ )o always preceded by an A^— particle equilibrium 
distribution function po- It includes the overlap function W{N), which vanishes for overlapping configurations (with 
at least one pair distance < cr), and W{N) = 1 for non-overlapping configurations. For the derivation of these and 
other properties we refer to the original literature (IS] . 

The only additional properties needed in this article are the conservation laws for the summational invariants = 
{l,^ri,vf}, where 

r±(zj)(a,+a,) = 0. (A.4) 

An important application of ljA.3|) in the present context is the interpretation of the Einstein-Helfand formulas 
(|II.8p with the time correlation function ((M(t) — M(0))^)^ for hard sphere fiuids. There is a caveat here. The 
phase function with A{t)B{t) and i > can only be expressed in hard sphere generators as e*'^+^(0)i3(0). It is not 
equal to (e'^+ A(0)(e*'^+i?(0)), because e*^+ is not a substitution operator, like the streaming operator e^^ for smooth 
interactions. One shows stationarity of an equal-time average as {A(t)B{t)) ^ — {e*^^ AB) ^ = {AB)^, where we have 
used (|A.3p and the relation (|A.3|I in the form e~*^-l = 1. 

With this in mind the Helfand formula (|II.8|I , say for the shear viscosity, can be expressed in hard sphere generators 



= UJ ^^-^^^ ^ -^(Afe'^+L+M)„. (A.5) 

Application of (jA.3|l and use of the identity, 

e-*^-=l-[ dse-'^-L_, (A.6) 







allows us to transform the Helfand formula for hard spheres into. 



Vv{t) = -^{ML+M)^ + ^ I ds{{L^M)e^^+L+M)^. (A.7) 







This relation is used in (IIII.2I) of the main text. 
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